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Rashba-split surface band is characterized by a one-to-one correspondence between the electron's 
momentum k and its spin orientation. Here we show that a similar correspondence between mo- 
mentum and orbital angular momentum (OAM) must exist on surface bands once the inversion 
symmetry is broken. The correspondence is valid even when there is no spin-orbit interaction. 
Tight-binding and first-principles calculations are presented to support our claim. As a method to 
detect such OAM-momentum correspondence, we propose the circular dichroism (CD) experiment 
using the angle-resolved photoemission (ARPES) setup. CD-ARPES experiment performed on Cu 
surface confirms the existence of chiral OAM. A new concept of "orbital Galvanic effect" is proposed. 
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Electrons in solids are classified by several quantum 
numbers, including their momentum, spin, and some- 
times orbital angular momentum (OAM). In some cases 
two or more degrees of freedom appear in coupled form 
as, for instance, in the Rashba phenomena where there 
exists a one-to-one correspondence between the electron's 
momentum and its spin[l]. Rashba effects are promi- 
nent in certain surface bands [2] or in a bulk insulator [3] 
lacking the inversion symmetry. In addition to inversion 
symmetry breaking (ISB), spin-orbit interaction (SOI) is 
the pre-requisite for the observability of Rashba-related 
phenomena. The relation of SOI to Rashba splitting has 
been discussed by several authors in the past [4], and re- 
cently by some of the present authors [5] . 

The implication of ISB on the Rashba splitting is an- 
ticipated on symmetry grounds. With the surface normal 
along the z-direction, one can write down a symmetry- 
allowed Hamiltonian Hji = Arz • (k x cr) in terms of 
the spin operator <t/2 and electronic momentum k. The 
Rashba energy scale Ar appearing in Hb, is on the or- 
der of the electrostatic potential barrier across the sur- 
face, multiplied by its width[5]. The chiral spin angular 
momentum (SAM) structure in momentum space follows 
as a direct consequence of the Rashba Hamiltonian Hb.. 
Such chiral SAM structure has been thoroughly docu- 
mented on surfaces of several metallic elements[6, 7] as 
well as on the surfaces of topological insulators [8] in re- 
cent years. Its potential as an effective source of spin 
current is also being actively investigated [7]. 

Upon closer inspection, however, one finds that other 
physical quantities with the same symmetry properties 
can take the place of cr in the Rashba Hamiltonian. It 
is conceivable, for instance, to replace cr by the orbital 



angular momentum operator L in a degenerate orbital 
system where L is an active degree of freedom. In the 
case of materials with strong spin-orbit interaction (SOI), 
the total angular momentum J ^ L-\- (l/2)cr may be the 
more appropriate choice [5]. Both scenarios would imply 
chiral OAM structure in nruch the same way that chiral 
SAM follows from the Rashba Hamiltonian. We show 
that chiral OAM in one-to-one correspondence with the 
electron's linear momentum is indeed a general conse- 
quence of ISB at the surface and should occur on surface 
bands even for materials with no SOI, placing OAM as 
the more generic feature of ISB than SAM is. 

The above symmetry consideration is further sup- 
ported by the analysis of a tight-binding (TB) model of 
two-dimensional monolayer of atoms. For simplicity we 
will consider an sp-orbital system with three degenerate 
p-orbitals forming bands. To emphasize the notion of 
ISB-induccd OAM better, we will first study the spinless 
case to show that OAM can arise without SOI and only 
later introduce SOI as a perturbation. The latter proce- 
dure is shown to recover the usual Rashba spin splitting. 
It will be shown that the phenomenon of chiral OAM is 
described by a two-dimensional massive Dirac Hamilto- 
nian. 

The tight-binding model of spinless p-orbitals can be 
constructed in terms of two Slater-Koster parameters Vi 
and V2 for a- and 7r-bonding amplitudes, respectively, 
and a third one, 7, representing the degree of ISB [4]. In 
real materials, 7 arises from the surface-normal electric 
field which breaks the inversion symmetry. Degeneracy 
of atomic p-orbital states is assumed in our model, which 
is justified in simple elements like Cu, Sb and Bi due to 
their weak crystal field splitting. We will write N for the 



number of sites in the lattice, and \i,p\) for the locahzed 
Wannier orbitals (A = x,y,z) at the atomic site r^. Then 
the tight-binding HamiUonian for the triangular lattice 
(essentially the same result obtains for square lattice) 
near the F-point (k = 0) in the momentum-space basis 
|k,pA) = N^^^'^J2i'^^^''^'\hP\) becomes 



akl + /3ky {a-l3)kxky 
H]^= I {a — f3)kxky aky + /3kl 
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Here, a = 3(3Fi - F2)/8, /3 = 3{Vi - 3t/2)/8 and P = 
k^. + ky. The lattice constant is taken to be unity. To 
diagonalize _ffk, it is convenient to choose a new set of 
basis vectors 



|I,k) = (fc,/fc)|p„k)-(fc,/fc)|p„k), 

|II,k) = {k,/k)\p,,k) + {ky/k)\Py,k), 

|III,k) = e-'^-\p,,-k), (2) 

k = |k|, e^'^^ — {kx +iky)/k. The state |I,k) remains 
decoupled at energy Eim = 3V^2 - 3Fi + 3(T4 - 3V2)k^/8, 
while I II, k) and |III, k) obey a reduced 2x2 Hamiltonian 



HoAM = -4/3/2x2 + M-'k" + -z • (7k X T - At'), (3) 



with A = Vi + 1^2 the bandwidth, M^^ == ( " „, . ,„ 

yu —31/2/2 

the effective mass tensor, and r the pseudo-spin matrix. 
Eigenstates in the leading order of 7/A are 



|2,k) 2. |II,k)- 
|3,k) ^ |III,k) 



»7(fca; - iky) 

2A 



III,k), 



il(kx + iky),„ , > 
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with energies E2M ^ HV2 - Vi) + 3((3Vi - T/2)/8 - 
72/4A)k2 and E^^y, ~ 6V2 + (37V4A - 3F2/2)k2. The 
0AM operator is given by the sum L = i^/N) J^i ^i 
where each Li acts on the Wannier state \i,p\) as the 
usual L — 1 angular momentum operator. The two 
bands obtained above carry nonzero 0AM as claimed 
(L+ = L"" + iLy): 



{2,k\L+\2,k)^'-^ikx+iky) 



-(3,k|L+|3,k). (5) 



Both bands possess a chiral pattern of 0AM whose 
strength is proportional to the ISB parameter 7 and are 
opposite between the two bands. 0AM for the decoupled 
state |I, k) is vanishing because it consists of Pa; andpy or- 
bitals only {Lz — ±1) and does not contain pz {Lz ~ 0). 
The reduced Hamiltonian (3) satisfies time-reversal in- 
variance. The size of 0AM gets reduced linearly with 



the momentum |k| and with the symmetry-breaking pa- 
rameter 7. By contrast, the magnitude of polarized spin 
in the Rashba model is independent of the Rashba pa- 
rameter Ar. 

Incorporating the spin degree of freedom at each 
atomic site as the Pauli matrix <Ti, we can examine the 
influence of the SOI interaction iJgo = (Aso/2) X^i ^i ' ""* 
as a perturbation for each band obtained above. The 
spin-orbit interaction energy Ago is assumed smaller than 
either the ISB energy 7 or the bandwidth A. Employ- 
ing the standard degenerate perturbation theory for the 
two spin states at a given k-vector, we find the resulting 
matrix elements {n,k, s\Hso\n,k, s') {s,s' =t,4.) form a 
2x2 Rashba-type Hamiltonian matrix 



± — Aso-r2; • k X (T 

2 A 



(6) 



for bands n = 2 and n == 3, respectively. There is no 
Rashba splitting for band 1 within the degenerate per- 
turbation theory. Combined with the analysis of chiral 
0AM, we obtain the following hierarchy of chiral angular 
momenta: First, there is chiral 0AM arising from ISB 
alone. At the next level, when a small SOI is present, 
each spin-degenerate OAM-carrying band is split into a 
pair of bands carrying opposite chiral SAMs. Both chiral 
structures are embodied in Dirac-like effective Hamiltoni- 
ans, with the energy scale for chiral 0AM dominant over 
that of chiral SAM by the factor A/A^o. Chiral 0AM, 
unlike SAM, can be substantial for surfaces consisting of 
light elements, and its strength controlled by applying an 
inversion-symmetry-breaking electric field externally. 

To ensure that chiral 0AM exists in a more realis- 
tic calculation, we performed first-principles local-density 
approximation (LDA) calculation for a Bi single layer 
forming a triangular lattice. The choice is inspired by 
Bi being a proto-typical p-orbital band material. An ex- 
ternal electric field of 3V/A perpendicular to the layer 
was imposed by hand to mimic the surface potential gra- 
dient without having the complication of dealing with 
the bulk states. The physically more relevant case of 
Bi layer is also considered [9], with LDA results easily 
adaptable to the single-layer case discussed below. To 
emphasize the relevance of ISB we chose to investigate 
the spin-degenerate by turning off SOI in the LDA cal- 
culation. The resulting electronic structure for spinless 
case consisting of three p-orbital-derived bands is shown 
in Fig. 1(a). As the external electric field is turned on, 
a level repulsion between the middle {E2 in Fig. 1(a)) 
and the bottom {E3 in Fig. 1(a)) band occurs as in- 
dicated by circles in Fig. 1(a). These two bands exhibit 
the chiral 0AM patterns with the maximum 0AM vector 
\{L)\ w 0.96fi as shown in Fig. 1(c) and 1(d), while the 
third one, shown in Fig. 1(b), carries much less 0AM 
around the F point. The 0AM chiralities of the two 
bands are opposite, in accordance with the previous TB 
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FIG. 1: OAM and CD from first-principles and tiglit-binding 
calculations of Bi monolayer without SOI. (a) LDA band 
structure for Bi monolayer of triangular lattice structure with- 
out SOI. Perpendicular electric field of 3V/A was imposed 
externally. Three dashed curves represent the tight-binding 
energy dispersions around the F point, (b)-(d) OAM vectors 
(green arrows) and CD signals (color backgrounds) for the 
three bands from the highest (-Bi) to lowest (-B3) energies over 
the whole Brillouin zone marked by a solid hexagon. Largest 
OAM has magnitude ~ 1ft for bands E2 and -E3. Incoming 
photon direction kph = +x is assumed for the CD calculation. 
(e)-(f) OAM obtained from the tight-binding calculation for 
two main OAM-carrying bands E2 and £3. 



analysis. An excellent fit of the LDA band structure 
near the T point was possible with the TB parameters 
Vi = --0.725 eV, V2 = -0.11 eV and 7 = 0.2623 eV 
(Fig. 1(a)). The good fit by the TB model is a testa- 
ment to the negligible crystal field effect in the elemental 
Bi. OAM values obtained by the TB analysis are also 
in excellent quantitative accord with the LDA results as 
can be seen by comparing Fig. l(c)-(d) (LDA) to Fig. 
l(e)-(f) (TB). The OAM magnitude is seen to decrease 
continuously upon approaching the F point in the LDA 
calculation (Fig. 1(c) and 1(d)) as predicted by the TB 
calculation, Eq. (5). 

Having established theoretically the existence of chi- 
ral OAM in inversion-asymmetric bands by a num- 



ber of methods, we turn to the question of its de- 
tection. Spin- and angle-resolved photoemission spec- 
troscopy (SARPES) has served to identify the chiral spin 
structure of the surface bands in the past [6, 7]. A simi- 
lar chiral structure for OAM as demonstrated here can- 
not, however, be detected by the same probe since chi- 
ral OAM exists even when SOI is very weak and spin 
degeneracy is nearly perfect. Circular dichroism (CD) 
refers to phenomena in which the physical response of 
a system to probing light depends systematically on the 
light polarization being left-circularly-polarized (LCP) or 
right-circularly-polarized (RCP). In ARPES experiment, 
in particular, incident lights of opposite helicities (RCP 
vs. LCP) might give rise to different scattering intensi- 
ties of photo-electrons. It can be shown, in fact, that such 
difference is a consequence of finite OAM polarization in 
k-space and can be used to detect its existence. 

One can formally define the CD- ARPES signal £)(k) 
as 



Z?(k) 



E.(/.^^^(k)-/.^^P(k)) 

E.(/«CP(k)+/LCP(k))' 



(7) 



Here /cr(k) refers to the probability of scattering from 
an initial occupied state at momentum k to a final 
photo-electron state of spin orientation a. Since we have 
spin-integrated ARPES in mind, the final-state spin is 
summed. Coupling to electromagnetic fields gives rise to 
a perturbed Hamiltonian Hi ^ p • A (p=momentum op- 
erator, A=vector potential). The transition amplitude 
from the initial state |/) to a final state |_F) is propor- 
tional to {F\Hi\I) ^ (_F|r ■ A|/), according to standard 
theory. Taking a plane- wave form for the final state wave 
function ippi^) = e'*'''"", where kp is the wave vector of 
the final-state photo-electron, and using the initial state 
wave function obtained from LDA, we can calculate the 
transition probability for RCP (A) and LCP (A*) incom- 
ing lights to get D{k). The in-plane component of k^ 
matches the momentum of band electrons (higher-order 
Umklapp processes are neglected) , while the normal com- 
ponent kp can be deduced from energy conservation. 

The incident photon direction fcph is chosen with its 
in-plane component along -\-x axis in calculating -D(k) 
for Bi monolayer band structure. The results are color- 
coded and superimposed on the OAM patterns in Fig. 
l(b)-(d). In the case of Bi monolayer the CD intensity 
proved largely insensitive to the choice of kp value, so 
we used kp — for the displayed result. One can make 
sense of the CD color patterns from a simple angular 
momentum conservation argument when k is close to the 
r point. Defining the projection of the electron angular 
momentum operator L onto the photon direction as /cph • 
i, incident RCP light would increase it by +fi. whereas 
LCP light reduces it by h. Therefore, loosely speaking, 
initial states with {L) oc —x {l=-l) is coupled to Z = 



(l = -2) channel of the final state by the RCP (LCP) 
light. Expansion of the final plane-wave state in terms 
of spherical harmonics produces a larger weight for / = 
than I = —2, thus I?(k) > for such an initial state. 
Following such consideration leads to the result £'(k) oc 
V-(k|i|k). 

In Fig. 1(c)- (d), one finds that such elementary con- 
sideration to be indeed in accord with the CD patterns 
near the F point. Far from it, though, the color pattern is 
no longer consistent with it and even gives out the oppo- 
site sign, implying that effects besides the conservation 
rule play a crucial role in determining the overall sign 
of £'(k). In principle, however, all such effects are cap- 
tured in the LDA calculation and should be compatible 
with ARPES experiment. We emphasize that CD signal 
vanishes where 0AM disappears, leaving little doubt that 
the measurement of CD is a direct indicator of local 0AM 
in momentum space. On the other hand, the existence 
of spin polarization in the band will not be captured by 
the proposed method even when it is prominent. 

Finally, we performed CD- ARPES experiment on Cu 
surface with a view to confirm the existence of chiral 
0AM on a realistic material with a small SOI. Because 
SOI is very small in Cu, it can be safely said that chi- 
ral 0AM structure is, without question, due to the ISB 
mechanism proposed here. The calculated band struc- 
ture of 30 Cu(lll) layers is shown in Fig. 2(a). Two 
surface-derived bands are identified as red curves. In the 
case of Cu calculation the inversion symmetry-breaking 
electric field is generated spontaneously through self- 
consistent electronic structure calculation, and not im- 
posed externally as in the Bi monolayer case. Despite 
these differences we found a clear, chiral 0AM pattern 
around the F point as shown in Fig. 2(b). Even though 
the make-up of the surface bands receive substantial con- 
tributions from both p- and d-orbitals, 0AM arises pre- 
dominantly from the d-orbital components. Using Eq. 
(7) we can also calculate r'(k) for Cu. Contrary to the 
model Bi calculation, £'(k) varies significantly with the 
choice of the final-state momentum kp. For Fig. 2(b) we 
used kp = 2.27A~^ consistent with the photo-electron 
energy in the actual experiment described in the next 
paragraph. The sign of -D(k) is opposite to the antic- 
ipation from the angular momentum conservation rule 
near the F point. It agrees, however, with the following 
dichroism experiment performed on Cu surface. 

The single, sharp experimental Cu(lll) surface band 
shown in Fig. 2(c) actually consists of two surface bands 
with opposite spins as found in LDA calculation, unre- 
solved because of the small energy spacing between them. 
We take the ARPES data twice, one with RCP and the 
second with LCP lights, and obtain £'(k) according to 
Eq. (7). The results, plotted in Fig. 2(d), cover the small 
occupied region in Brillouin zone with the Fermi mo- 
mentum ss 0.2A~^. The CD-ARPES experiment shows 
pronounced CD behavior that reaches a maximum value 
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FIG. 2: OAM and CD of Cu surface band, (a) First-principles 
electronic band structure of Cu(lll) 30-layer slab. Red 
dashed lines indicate the two surface bands, (b) Calculated 
OAM and CD (assuming the in-plane incident photon direc- 
tion kph oc +x) corresponding to outer energy surface bands. 
Inner surface band shows basically the same OAM and CD 
patterns. Maximum OAM vector is ~ Q.07fi. (c) ARPES- 
measured surface band of Cu along ky = 0. (d) Measured 
CD-ARPES result D(k). Incident photon direction has in- 
plane component along +x. The momentum window in Fig. 
(b) matches that of Fig. (d) . 



of « 20% at the Fermi surface, where the OAM is also 
expected to be the largest. Most importantly, the sign 
of CD matches that in the LDA calculation, which is a 
highly non-trivial result given the reversal of the CD sign 
discussed earlier. The calculated D(k) varies between - 
0.3 and -1-0.3, somewhat greater than the 20% variation 
found experimentally. Thermal effects in the experiment 
may have contributed to the smearing of CD signal. 

Our prediction, and its observation on Cu surface, of 
chiral OAM is a direct analogue of the well-known chiral 
SAM from Rashba-split surface bands. Likewise, much 
of the physical consequences of chiral SAM is expected 
to have orbital analogues. For instance, spin polariza- 
tion can drive electrical current through spin-momentum 
locking in a phenomenon known as the spin Galvanic 
effect [10]. Its inverse, the spontaneous spin polarization 
in the current-carrying state of Rashba-split bands, was 
proposed theoretically by Edelstein[ll] and confirmed ex- 
perimentally in recent years[12]. Having identified OAM 
as the more fundamental physical quantity associated 
with ISB, we propose that both these spin-related effects 
must have orbital analogues. 
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